Linear and nonlinear optical characteristics of the Falicov-Kimball model 
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We calculate the linear and nonlinear optical properties of the Falicov-Kimball model for a mixed- 
valent system within the self-consistent mean-field approximation. Second-harmonic generation can 
only occur if the mixed-valent state has a buih-in coherence between the itinerant d-electrons and 
the locahzed /-holes. By contrast, second-harmonic generation cannot occur for solutions of the 
model with /-site occupation as a good quantum number. As an experimental test of coherence in 
mixed-valent compounds we propose a measurement of the dynamic second-order susceptibility. 



The Falicov-Kimball 0| (FK) model has been used ex- 
tensively for the mixed-valent compounds, heavy fermion 
systems, and associated metal-insulator transitions. The 
FK model accounts for a band of itinerant d-electrons 
and localized /-orbitals and intrasite Coulomb interac- 
tion between the d- and /-electrons. A d-f hybridization 
term may or may not be added to the model. The the- 
oretical solutions for the ground state of the FK model 
can be divided into two classes. On the one hand, solu- 
tions that treat the occupation of an /-electron on a site 
as a good quantum number do not have a built-in 
coherence between c?-electrons and /-holes. On the other 
hand, solutions such as the self-consistent mean-field so- 
lution Q and the electronic polaron ||] do have a built-in 
coherence between c?-electrons and /-holes. 

We report here the nonlinear optical responses of these 
two classes of solutions. Solutions of the model with a 
built-in coherence can sustain second-harmonic genera- 
tion. Solutions with classical /-electron site distributions 
cannot. Therefore, we propose the measurement of the 
second-order susceptibility of a mixed-valent compound 
as a test to distinguish between these theories. The exis- 
tence of such second-harmonic generation due to coher- 
ence in the ground state would, of course, be of interest 
in its own right as a manifestation of strong electron cor- 
relation. 

Four- Wave-Mixing (FWM) spectroscopy has become 
a powerful tool for studying coherence in semiconduc- 
tor systems 0. In a three-beam FWM experiment, two 
incoming beams of wavevectors ki and k2 set up a tran- 
sient polarization grating. The third incoming beam of 
wavevector ka diffracts off the grating to produce the 
outgoing signal in the direction k4 = ks + k2 — ki . Be- 
ing a third-order process, FWM is allowed in media with 
or without inversion symmetry. We pose the question: 
what happens if the state being probed already has a 
polarization built into it instead of being created artifi- 
cially by optical pumping? An example of such a system 
is the self-consistent mean-field (SCMF) solution of the 
FK model resulting in the Bose-Einstein condensation of 
d-f excitons. 

As shown below, the built-in polarization leads to a 
nonlinear optical response to second order in the exter- 
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nal field. The built-in polarization replaces one of the 
external fields of the three-beam FWM experiment. The 
mixed-valent system has a nonvanishing susceptibility 
X^^^(2a;, cj, cj) for second-harmonic generation. In crys- 
tals with inversion symmetry, second-harmonic genera- 
tion is forbidden under the electric-dipole approxima- 
tion. In the mixed-valent system, the built-in polar- 
ization breaks the inversion symmetry, allowing second- 
harmonic generation to take place. The built-in polar- 
ization also means that the system is ferroelectric. Our 
calculation of the concomitant dielectric behavior will be 
reported elsewhere. There are reports of usually large 
dielectric constants in mixed valent semiconductors [Q. 
Because of the problem of residual carriers in dielectric 
measurements, we feel that the nonlinear optical effect 
might be a clearer test. 

We present a calculation of the linear and second- 
harmonic susceptibilities of a model mixed-valent sys- 
tem within the SCMF approximation. The magnitude 
of the second-harmonic output signal is directly propor- 
tional to the built-in coherence A. The Coulomb interac- 
tion between the optically excited quasiparticles greatly 
enhances the second-harmonic conversion efficiency at 
a; = A (one half the energy gap 2A). 

Ignoring the electron spin, the FK Hamiltonian is 

^ - E ^k^'^k + 4 E /k/k + E(^k4/k + h.C.) 
k k k 

+^ E 4+crfk/^,_,/k'. (1) 

k,k',q 

Here dj^ (/^) creates a d-(/-)electron of momentum k and 
energy Ek (£'/)■ The parameters U and Vk are the direct 
interaction and the hybridization between the d- and /- 
electrons, and N is the number of sites. We consider a 
model system with a d-band and /-level arising from d- 
and /-orbitals on the same site. The d-band has band- 
width 2W and a constant density of states po = 1/(2M^). 

The SCMF solution is analogous to the BCS theory 
of superconductivity except that the pairing now occurs 
between a c?-electron of momentum k and a f-hole of 
momentum — k. The ground state is {ip) = nk(^k + 
Vkdl^fk)\0), where |0) is the state with no /-holes (the 
normal state), and (uk) is the probability amplitude 
for the pair state (k, — k) to be occupied (unoccupied). 

A key feature of \tp) is that it is a state of broken 
inversion symmetry. If the crystal is invariant under in- 
version with respect to a d-f site, the inversion symme- 
try is broken by the pairing of d-states of even parity 
with /-states of odd parity. Applying the inversion J on 

yields the state J\tp) = nk('^k ~ '''k^k/k)|0): which 
is orthogonal to lip) and has the same energy except in 
the case U — when the two states are the same. The 
degenerate states and J\4>) have built-in polariza- 
tions of opposite directions, for the polarization operator 
P = (X^k/^^k/k + h.c.)/ri where O is the system vol- 
ume. We take the interband dipole matrix element fi to 
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be independent of k. The correct ground state is selected 
by lifting the degeneracy with an infinitesimal external 
electric field E, and choosing the lower energy state. The 
consequent breaking of the inversion symmetry is what 
allows second-harmonic generation to take place. 

The built-in polarization defines a direction in space, 
which we call the z-axis. (Without crystal-field terms, 
the z-axis has no definite orientation with respect to 
the crystal axes.) Since fi^ is real, Pi^"* — Nfi^i^ + 
A*) /{^1U), where A is the built-in coherence. The built- 
in coherence is determined self-consistently from 



where 2Ei, = ^(ck - Ef)^ + 4|A - Vu]^ is the quasi- 
electron-hole pair excitation energy g. Eq. (g) is 
Eq. (11) of Ref. [Q with a k-dependent hybridization. 
If the crystal is invariant under inversion, the hybridiza- 
tion must satisfy = — Vk. Then, since Vk is purely 
imaginary and odd in k, the imaginary part of A van- 
ishes due to the cancellation of terms with ±k. The real 
part of A is given by the BCS gap equation 



with 2Ei, = ^(ck - Ef)^ + 4A2 + 4|ykP- Calculation 
shows that sufficiently strong Vk can destroy the gap. In 
the following we consider the limit where Vk is neligible. 

A is the order parameter of the valence transition. 
When the /-level is far below the d-band, the system 
is in the normal state (A = 0). As the /-level is moved 
upward past a critical value (in a real material this is 
done by applying pressure or alloying), the system un- 
dergoes a transition into the mixed- valent state (A > 0). 
In the mixed- valent state, the /-level occupancy nj lies 
between zero and one. A reaches a maximum at the half- 
filling point Ef = 0. Electron- hole symmetry requires 
A{-Ef) = A{Ef) and nf{-Ef) = 1 - nf{Ef). 

We first consider the linear absorption spectrum of the 
mixed- valent system. The SCMF solution predicts an en- 
ergy gap 2 A. Far-infrared optical measurements [p|-pT|, 
as well as electron tunneling experiments jl2j, show an 
energy gap of several meV in a number of mixed-valent 
compounds. The crucial difference between the super- 
conductor and the mixed-valent system is this: in the 
superconductor, the pairing occurs between two elec- 
trons, whereas in the mixed-valent system, the pairing 
occurs between an electron and a hole. This has im- 
portant consequences for the coherence factors that en- 
ter the response of both systems to different external 
probes. For example, the coherence factor entering the 
optical absorption of the mixed-valent system is the same 
as the coherence factor entering the acoustic attenuation 
of the superconductor. The interaction of the mixed- 
valent system with the electromagnetic field is treated in 
the electric-dipole approximation. Only the component 
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of the electric field along the symmetry-breaking z-axis 
couples to the channel in which the pairing takes place. 
We ignore the response of the remaining optical chan- 
nels. Second-harmonic generation can only occur in the 
symmetry-breaking channel. 

We have calculated the linear susceptibility Xz'z both 
from the Kubo formula and from the optical Bloch equa- 
tions. The pseudo-spin picture gives a nice physical de- 
scription of the linear and nonlinear responses as pre- 
cessional modes of the pseudo-spin Sk. For a given k, 
the pseudo-spin corresponds to the dk- and the /k-state 
as a two-level system. The equations of motion for the 
pseudo-spin are the optical Bloch equations 

Sk = (Hk - Mk) X Sk. (4) 

Here Hk — (— 2^2i?2, 0, ek ^ ^/) is the external "mag- 
netic" field, and Mk — jj J2k is the pseudo-magneti- 
zation. The symbol x means the vector cross product. 

To calculate the linear susceptibility we expand Sk, 
Hk and Mk to first order in the electric field E^. From 
Eq. (^) to zeroth order we obtain an equation for the 

stationary pseudo-spin s'^^ . The built-in coherence tilts 

S^""* away from the negative z-axis. The tilting angle 

is ^k = arccos(wk ~ ^^k)- F'or A real, S^'^^' lies in the x-z 
plane. The electric field causes the pseudo-spin to precess 
around the stationary direction. With the precession axis 
tilted away from the z-axis, the field causes variations in 
all three Cartesian components of Sk. A simpler descrip- 
tion is obtained in the spherical polar coordinate system. 
The stationary direction is the radial unit vector e^. The 
precession is decomposed into components along the po- 
lar and azimuthal unit vectors and e^. The equations 

of motion for the polar and azimuthal components Sg^l^ 
and 5*^^^ are 

si]l~'2E^sl'l + Ml'l = F^% (5) 
si]l + 2E^sl,]l-Ml,'l^Fi% (6) 

The driving terms are Fg^^ ~ and F^^\^ = 2fizEz cos6'k. 
The linear susceptibility is xi^z — Pz^^ /E^, where Pi^"* = 
(-^/^z X^k "^e^k '^'^s6'k)/i^ is the linear polarization. 

For a separable interaction potential, Eqs. (||) and (||) 
can be solved analytically. The linear susceptibility is 

.(1) . ^ ( m 1^ 



where 



"^^"^^ 2E^{uj~2Ej,){uj + 2Ei,)' 

, , U ^— V £k — Ef 
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For the simple model system, A{u!) and B{uj) can be ex- 
pressed in terms of elementary functions. The poles of 
Xzz give the collective excitation energies of the mixed- 
valent system. The denominator in Eq. (0) vanishes when 
(w^ - 4A2)A2(a;) - B'^{lo) = 0. The zero-frequency Gold- 
stone mode is a consequence of the arbitrariness of the 
phase of A (in the absence of hybridization). In the 
pseudo-spin picture, the Goldstone mode corresponds 
to rotating all pseudo-spins around the z-axis over the 
same angle (j). Since this does not change the total en- 
ergy, the Goldstone mode has zero frequency. In Eq. (0), 
the Goldstone mode does not appear to contribute to 
the linear optical response since the pole at a; = is 
canceled by a factor of uj in the numerator. There are 
no exciton-like collective modes within the energy gap. 
When u! < 2A, the functions A{uj) and B{uj) are purely 
real, so {uj^ - 4A'^)A'^{uj) - B'^{uj) < 0. 

The absorption spectrum is given by the imaginary 
part of xiV- When the /-level lies inside the d-band 
{\Ef \ < W), the absorption spectrum has a threshold 
singularity at w = 2A. When \Ef\ < W, the singular- 
ity is e~^/^6'(e), and when \Ef \ = W the singularity is 
e"i/2ln"^(e)6'(e), where e = lu - 2A. When the /-level 
lies outside the c?-band the singularity is cut off because 
the energy gap is larger than 2A. The singularity is due 
to the final-state Coulomb interaction between the op- 
tically excited quasiparticles. In the single-quasiparticle 
picture, the absorption spectrum rises continuously from 
zero according to e^/^0(e). The singularity is not an ar- 
tifact of the simple model, and should be observable in 
real materials. 

(2) 

We calculate the second-harmonic susceptibility Xzzz 
from the optical Bloch equations by expanding the 
pseudo-spin and the pseudo-magnetization to second or- 
der in the perturbing electric field Ez- The equations of 

(2) (2) 

motion for the second-order components Sg I. and 
have the same form as Eqs. and (|6|), except with more 
complicated driving terms. A very important observa- 
tion is that all driving terms are directly proportional 
to A. When A = 0, the second-harmonic susceptibility 
vanishes identically. 

In addition, the second-order fluctuations have a non- 

(2) 

zero radial component S^. ^ . The motion is no longer a 
regular precession: the pseudo-spin nutates during the 
precession. (Nutation is the up-and-down motion of the 
precession axis.) The nutation frequency is twice the pre- 
cession frequency. The second-harmonic susceptibility is 

= pP/EI where = [iV/i, ^^(5^'^ cos 0k + 

Sl.^sm0i^)]/Q is the second-order polarization. 

For a separable interaction potential, an analytic so- 
lution for the second-harmonic response is possible in 
principle. However, the large number of driving terms 
presents a considerable challenge. We have instead ap- 
proached the problem numerically. This is done in anal- 
ogy with the classical mechanics treatment of forced os- 
cillations. One first solves for the motion in normal co- 
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ordinates, and then takes linear combinations to obtain 
the motion in the original coordinates. The results of 
the calculation are shown in Fig. pT. The figure shows the 

(2) 

amplitude \xzzz{2uj,uj,uj)\ of the second-harmonic sus- 
ceptibility as a function of the photon energy lo, for 
various values of Ef. The important features are: (1) 
The second-harmonic amplitude is directly proportional 
to the amount of coherence A built into the mixed- valent 
system. (2) When the /-level lies inside the d-band, 
the second-harmonic conversion efficiency is strongly en- 
hanced at a; — A, and less strongly at u; = 2 A. The 
first feature shows that second-harmonic generation can 
be used as a test of the validity of the SCMF solution in 
real mixed-valent materials. The second feature distin- 
guishes the single-quasiparticle treatment of the second- 
harmonic response from the self-consistent mean-field 
treatment. Like the threshold singularity in the case of 
linear reponse, the enhancement of the second-harmonic 
conversion efficiency is due to the final-state Coulomb 
interaction between the optically excited quasiparticles. 

As an experimental test of coherence in mixed-valent 
compounds we propose a measurement of their second- 
harmonic susceptibility x*^^) (2w, oj, oj). Consider for ex- 
ample SmBg. The crystal structure of SmBg has cubic 
symmetry, with Be octahedra at the body center, and Sm 
ions at the corners of a conventional bcc unit cell with lat- 
tice constant a = 4. 13 A. The crystal has inversion sym- 
metry at the bcc lattice points. Through measurements 
of the ionic radius, the valence of the Sm ion in SmBg is 
found to be 2.53, almost halfway between 2 and 3, so that 
the /-level lies near the center of the conduction band. 
The measured far-infrared absorption spectrum P-[Tl]] of 
SmBg can be interpreted in accordance with the SCMF 
solution. In Figure ^ we compare the mean-field and 
single-quasiparticle results for the linear susceptibility to 
experimental data on SmBg taken from Ref. The 
data show an energy gap around 2 A = 4meV, and a 
sharp peak at threshold. The mean-field theory fits the 
data very well in the threshold region, whereas the single- 
quasiparticle theory gives a qualitatively wrong threshold 
behaviour. Away from threshold, discrepancies between 
mean-field theory and experiment occur because of our 
simple model density of states. Further experimental in- 
dication for the validity of the SCMF solution in SmBg is 
provided by the electron tunneling spectrum iQ, which 
can be explained by analogy with Giaever tunneling in 
superconductors. 

The second-harmonic generation for the other types of 
solutions of the FK model is now discussed briefly. If 
the /-occupancy at each site is a good quantum number 
(0 or 1), any solution |^,^, homogeneous or inhomoge- 
neous, will not give rise to second-harmonic generation. 
For the electronic polaron solution the nonzero co- 
herence yields a second-harmonic generation. We found 
A = 3.28W for the parameter values given in Ref. [^. 
To distinguish between the exciton condensation solu- 
tion and the electronic polaron requires an investigation 
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of the quantitative difference of their linear and nonlinear 
optical spectra. Such a theoretical study will be left for 
the future. 

In conclusion, we calculated the linear and nonlin- 
ear optical characteristics of the Falicov-Kimball model 
within the SCMF approximation. We found that the 
second-harmonic susceptibility is directly proportional to 
the amount of coherence A built into the mixed-valent 
system. We also found that the final-state Coulomb inter- 
action leads to a threshold singularity in the absorption 
spectrum, and strongly enhances the second-harmonic 
conversion efficiency at a; = A. As an experimental test 
of the validity of the SCMF solution in real mixed-valent 
materials we propose a measurement of x^'^'' {"^^ t ^ t ^) ■ 
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FIG. 1. Amplitude |xLz(2c<;, w,c<;)| of the second-har- 
monic susceptibility as a function of the photon energy w, 
for various values of the /-level Ef. The dash-dotted Une 
shows the phase of x^zzz{'2uJ,uJ,Ld) for Ej = — l.OVK. The 
Coulomb repulsion is ?7 = 3.0 W and the hybridization is 
Vk = 0. The ampUtude is given in units of iV/i^/(2nW^). 
For the parameter values given for the soUd Une in Fig. 2, 
Nnll{2Q.W'^) = 82nmV-^ 

FIG. 2. Comparison of the mean-field (solid line) and sin- 
gle-quasiparticlc (dash-dottod line) results for the imaginary 
part of the linear susceptibility x^zz [fjj) of SmBe to experi- 
mental data taken from Ref. 10 (diamonds). The /-level is 
Ef = 0, the bandwidth is W = 40meV, the Coulomb re- 
pulsion \s U = 0.38 W , and the hybridization is Vk = 0. The 
photon energy uj was given a small imaginary part 5 = Q.QIW . 
The interband dipole matrix element is jiz = ^Aeao for the 
solid line, and = 5.0 eao for the dash-dotted line. 
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